function [skewness,kurtosis] = skewnessAndKurtosis(data)

T  = size(data,1);
mu = mean(data,1);
y  = data - repmat(mu,T,1);
y2 = mean(y.^2,1);
y3 = mean(y.^3,1);
y4 = mean(y.^4,1);

skewness = y3./y2.^1.5;
kurtosis = y4./y2.^2;


end

